Chapter 7 



Hamiltonian dynamics 



Conservative mechanical systems have equations of mo- 
tion that are symplectic and can be expressed in Hamiho- 
nian form. The generic properties within the class of sym- 
plectic vector fields are quite different from those within 
the class of all smooth vector fields: the system always 
has a first integral ("energy") and a preserved volume, and 
equilibrium points can never be asymptotically stable in 
their energy level. 



— John Guckenheimer 




ou MIGHT THINK that the strangeness of contracting flows, flows such as the 



Rdssler flow of figure 2.6 is of concern only to chemists or biomedical 



engineers or the weathermen; physicists do Hamiltonian dynamics, right? 
Now, that's full of chaos, too! While it is easier to visualize aperiodic dynam- 
ics when a flow is contracting onto a lower-dimensional attracting set, there are 
plenty examples of chaotic flows that do preserve the full symplectic invariance of 
Hamiltonian dynamics. The whole story started with Poincare's restricted 3-body 
problem, a realization that chaos rules also in general (non-Hamiltonian) flows 
came much later. 

Here we briefly review parts of classical dynamics that we will need later 
on; symplectic invariance, canonical transformations, and stability of Hamiltonian 
flows. If your eventual destination are applications such as chaos in quantum 
and/or semiconductor systems, read this chapter. If you work in neuroscience 
or fluid dynamics, skip this chapter, continue reading with the billiard dynamics 
of chapter 8 which requires no incantations of symplectic pairs or loxodromic 
quaitets. 
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7.1 Hamiltonian flows 

(P. Cvitanovic and L.V. Vela-Arevalo) 

An important class of flows are Hamiltonian flows, given by a Hamiltonian appendix C 
H{q, p) together with the Hamilton's equations of motion remark 2.1 

dH dH 

opi dqi 

with the d = 2D phase-space coordinates x split into the configuration space 
coordinates and the conjugate momenta of a Hamiltonian system with D degrees 
of freedom (dof): 

■^ = (q,p), (1 = (qi,q2,---,qD), p = (pi,P2,---,pd)- (7.2) 

The equations of motion (7.1) for a time-independent, D-dof Hamiltonian can be 
written compactly as 

Xi = coijHjix) , Hj{x) = ^H{x) , (7.3) 

OXj 

where = (q,p) e Al is a phase-space point, and the a derivative of (•) with 
respect to Xj is denoted by comma-index notation (•),;, 



(^ = \\ \ \, (7.4) 



is an antisymmetric \dxd\ matrix, and I is the [DxD] unit matrix. 

The energy, or the value of the time-independent Hamiltonian function at the 
state space point x = (q, p) is constant along the trajectory x{t), 



d dH dH 

-H{q{t),p{t)) = —qi{t) + ^PAt) 
dt oqi a Pi 

dHdH dHdH ^ 
dqi dpi dpi dqi 



so the trajectories lie on surfaces of constant energy, or level sets of the Hamilto- 
nian {(q,p) : H{q,p) - E]. For 1-dof Hamiltonian systems this is basically the 
whole story. 
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Figure 7.2: A typical coUinear helium trajectory in 
the [ri, plane; the trajectory enters along the ri-axis 
and then, like almost every other trajectory, after a few 
bounces escapes to infinity, in this case along the ri- 
axis. In this example the energy is set to /f = £ = -1, 
and the trajectory is bounded by the kinetic energy = 
line. 




Example 7.1 Unforced undamped Duffing oscillator: When the damping term 
is removed from the Duffing oscillator (2.8), the system can be written in l-lamiltonian 
form, 




(7.6) 



This is a 1 -dof l^amiltonian system, with a 2-dimensional state space, the plane (q, p). 
The l-lamilton's equations (7.1) are 

q = p, p = q-q^ . (7.7) 

For 1-dof systems, the 'surfaces' of constant energy (7.5) are curves that foliate the 
phase plane (q, p), and the dynamics is very simple: the curves of constant energy are 
the trajectories, as shown in figure 7.1. 



Thus all 1-dof systems are integrable, in the sense that the entire phase plane 
is foliated by curves of constant energy, either periodic, as is the case for the 
harmonic oscillator (a 'bound state'), or open (a 'scattering trajectory'). Add one example B.1 
more degree of freedom, and chaos breaks loose. 



Example 7.2 CoUinear helium: In the quantum chaos part of ChaosBook . org we 
shall apply the periodic orbit theory to the quantization of helium. In particular, we will 
study coUinear helium, a doubly charged nucleus with two electrons arranged on a line, 
an electron on each side of the nucleus. The l-lamiltonian for this system is 

1,1,22 1 

H-:,pI+^pI + —— ■ (7.8) 

2 2 n r2 ri+ r2 
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Collinear helium has 2 dof, and thus a 4-dimensional phase space M, which energy 
conservation foiiates by 3-dimensionai constant energy hypersurfaces. In order to vi- 
sualize it, we often project the dynamics onto the 2-dimensional configuration plane, 
the {ri,r2), ri > quadrant, figure 7.2. It looks messy, and, indeed, it will turn out to 
be no less chaotic than a pinball bouncing between three disks. As always, a Poincare 
section will be more informative than this rather arbitrary projection of the flow. The 
difference is that in such projection we see the flow from an arbitrary perspective, with 
trajectories crisscrossing, in a Poincare section the flow is decomposed into intrinsic 
coordinates, a pair along the marginal stability time and energy directions, and the rest 
transverse, revealing the phase-space structure of the flow. 

Note an important property of Hamiltonian flows: if the Hamilton equations 
(7.1) are rewritten in the 2D phase-space form i,- - v/(jc), the divergence of the 
velocity field v vanishes, namely the flow is incompressible, V-v = (9,v,- = ojiHjj = 
0. The symplectic invariance requirements are actually more stringent than just 
the phase-space volume conservation, as we shall see in sect. 7.3. 

Throughout ChaosBook we reserve the term 'phase space' to Hamiltonian 
flows. A 'state space' is the stage on which any flow takes place. 'Phase space' 
is a special but important case, a state space with symplectic structure, preserved 
by the flow. For us the distinction is necessary, as ChaosBook covers dissipative, 
mechanical, stochastic and quantum systems, all as one happy family. 



7.2 Symplectic group 

Either you're used to this stuff... or you have to get used 
to it. 

— Maciej Zworski 
A matrix transformation g is called symplectic, 

g'cog = co, (7.9) 

if it preserves the symplectic bilinear form {x\x) = x^tox, where denotes the 
transpose of g, and co is a. non-singular [2D x 2D] antisymmetric matrix which 
satisfies remark 7.3 

co'^^-oj, a? = -l. (7.10) 

While these are defining requirements for any symplectic bilinear form, oj is often 
conventionally taken to be of form (7.4). 

Example 7.3 Symplectic form for D - 2: For two degrees of freedom the phase 
space is 4-dimensional, x - {q\,qi,p\,pi) , and the symplectic 2-form is 



f 





1 


^ 











1 


-1 











I 


-1 





J 



(7.11) 
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The symplectic bilinear form <x'' V'^^) is tlie sum over tine areas of tine paraiieiepipeds 
spanned pairwise by components of tiie two vectors, 

iW) = = - q?p['') + (q^pf - ifPi) ■ (7-12) 

it is this sum over oriented areas (not the Euclidean distance between the two vectors, 
|j,(2) _ f^gf ,-g preserved by the symplectic transformations. 



If g is symplectic, so is its inverse , and if g\ and gj are symplectic, so 
is their product g2g\- Symplectic matrices form a Lie group called the symplec- 
tic group Sp((i). Use of the symplectic group necessitates a few remarks about 
Lie groups in general, a topic that we study in more depth in chapter 10. A Lie 
group is a group whose elements ^(0) depend smoothly on a finite number A'^ of 
parameters (pa- In calculations one has to write these matrices in a specific basis, 
and for infinitesimal transformations they take form (repeated indices are summed 
throughout this chapter, and the dot product refers to a sum over Lie algebra gen- 
erators): 



g(6<p)^l+5(t}-T , 50 eR^, 1501 «1, (7.13) 



where {Ti,T2 • • • ,TAr), the generators of infinitesimal transformations, are a set 
of A'^ linearly independent \dxd^ matrices which act linearly on the (i-dimensional 
phase space M. The infinitesimal statement of symplectic invariance follows by 
substituting (7.13) into (7.9) and keeping the terms linear in 6(p, 

Tjco + (JTa = 0. (7.14) 

This is the defining property for infinitesimal generators of symplectic transfor- 
mations. Matrices that satisfy (7. 14) are sometimes called Hamiltonian matrices. 
A linear combination of Hamiltonian matrices is a Hamiltonian matrix, so Hamil- 
tonian matrices form a linear vector space, the symplectic Lie algebra sp{d). By 
the antisymmetry of oj, 

(cuT,y = wT„. (7.15) 

is a symmetric matrix. Its number of independent elements gives the dimen- 
sion (the number of independent continuous parameters) of the symplectic group 
Sp(fif), 

N = d{d+\)l2 = D{2D+l). (7.16) 

The lowest-dimensional symplectic group Sp(2), of dimension A'^ = 3, is isomor- 
phic to SU(2) and S0(3). The first interesting case is Sp(3) whose dimension is 
= 10. 
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It is easily checked that the exponential of a Hamiltonian matrix 



g-e^T (7.17) 



is a symplectic matrix; Lie group elements are related to the Lie algebra elements 
by exponentiation. 



7.3 Stability of Hamiltonian flows 



Hamiltonian flows offer an illustration of the ways in which an invariance of equa- 
tions of motion can affect the dynamics. In the case at hand, the symplectic in- 
variance will reduce the number of independent Floquet multipliers by a factor of 
2 or 4. 



7.3.1 Canonical transformations 

The evolution of /' (4.5) is determined by the stability matrix A, (4.9): 

j^fix) - A{x)J'{x) , Aijix) - coik Hkjix) , (7.18) 

where the symmetric matrix of second derivatives of the Hamiltonian, H^^n = 
dkdnH, is called the Hessian matrix. From (7.18) and the symmetry of H^kn it 
follows that for Hamiltonian flows (7.3) 

A'^ (jj + ojA = . (7.19) 

This is the defining property (7.14) for infinitesimal generators of symplectic (or 
canonical) transformations. 

Consider now a smooth nonlinear coordinate change form j,- = hi{x) (see 
sect. 2.3 for a discussion), and define a 'Kamiltonian' function K{x) - H{h{x)). 
Under which conditions does K generate a Hamiltonian flow? In what follows we 
will use the notation d~j = d/dyj, sij = dhi/dxj. By employing the chain rule we 
have that 

Kj - Hjsjj (7.20) 

(Here, as elsewhere in this book, a repeated index implies summation.) By virtue 
of (7.1), diH = -oJimym, so that, again by employing the chain rule, we obtain 

^ij^j^ - ~^ijSi,l^lmSm,nXn (7-21) 
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comple: 



Figure 7.3: Stability exponents of a Hamiltonian equi- 
librium point, 2-dof. 



degenerate saddle 



generic center 



saddle-center 



real saddle 

i \ (2) 



degenerate center 



The right hand side simplifies to x, (yielding Hamiltonian structure) only if 



'O^ijSljOJlynSfyin — 6 in 



{121) 



or, in compact notation, 



(7.23) 



which is equivalent to the requirement (7.9) that dh is symplectic. h is then called 

a canonical transformation. We care about canonical transformations for two 

reasons. First (and this is a dark art), if the canonical transformation h is very example B.1 

cleverly chosen, the flow in new coordinates might be considerably simpler than 

the original flow. Second, Hamiltonian flows themselves are a prime example of 

canonical transformations. 

Dream student Henriette Roux: "I hate these s„, „. Can't you use a more sensible 
notation?" A: "Be my guest." 



Example 7.4 Hamiltonian flows are canonical: For Hamiltonian flows it follows 
from (7.19) that ^ (J^coJ) = 0, and since at the initial time J^ixo) - 1, Jacobian matrix 
is a symplectic transformation (7.9). This equality is valid for all times, so a l-lamilto- 
nian flow f'(x) is a canonical transformation, with the linearization dyf'(x) a symplectic 
transformation (7.9): For notational brevity here we have suppressed the dependence 
on time and the initial point, J - J'{xo). By elementary properties of determinants it fol- 
lows from (7.9) that l-lamiltonian flows are phase-space volume preserving, \detJ\ = 1 . 
The initial condition (4.9) for J is = 1, so one always has 



dety = +1 



(7.24) 



7.3.2 Stability of equilibria of Hamiltonian flows 



For an equilibrium point Xq the stability matrix A is constant. Its eigenvalues 
describe the linear stability of the equilibrium point. A is the matrix (7. 19) with 
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real matrix elements, so its eigenvalues (the Floquet exponents of (5.1)) are either 
real or come in complex pairs. In the case of Hamiltonian flows, it follows from 
(7.19) that the characteristic polynomial of A for an equilibrium Xq satisfies 

det (A - Al) - det {aj~\A - Al)aj) - det {-ojAco - Al) 

= det (A"' + Al) - det (A + Al) . (7.25) 

That is, the symplectic invariance implies in addition that if A is an eigenvalue, 
then -A, A* and -A* are also eigenvalues. Distinct symmetry classes of the Floquet 
exponents of an equilibrium point in a 2-dof system are displayed in figure 7.3. 
It is worth noting that while the linear- stability of equilibria in a Hamiltonian 
system always respects this symmetry, the nonlinear stability can be completely 
different. 



7.4 Symplectic maps 

So far we have considered only the continuous time Hamiltonian flows. As dis- 
cussed in sect. 4.4 for finite time evolution mappings, and in sect. 4.5 the iterated 
discrete time mappings, the stability of maps is characterized by eigenvalues of 
their Jacobian matrices, or 'multipliers.' A multiplier A = A(;t:o, t) associated to 
a trajectory is an eigenvalue of the Jacobian matrix /. As / is symplectic, (7.9) 
implies that 

-cor CO, (7.26) 

so the characteristic polynomial is reflexive, namely it satisfies 

det {J - Al) - det (/"^ - Al) - det {-coj'^a) - Al) 

- det(7~'-Al) - det(7"')det(l- A7) 
= A2^det(/- A-4). (7.27) 

Hence if A is an eigenvalue of /, so are 1 / A, A* and 1 / A* . Real eigenvalues 
always come paired as A, 1 /A. The Liouville conservation of phase-space vol- 
umes (7.24) is an immediate consequence of this pairing up of eigenvalues. The 
complex eigenvalues come in pairs A, A*, |A| = 1, or in loxodromic quartets A, 
1/A, A* and 1/A*. These possibilities are illustrated in figure 7.4. 

Example 7.5 Hamiltonian Henon map, reversibility: By (4.44) the Henon 

map (3.17) for b - -I value is the simplest 2-dimensional onentation preserving area- 
preserving map, often studied to better understand topology and symmetries of Poincare 
sections of 2 dof Hamiltonian flows. We find it convenient to multiply (3.18) by a and 



section 5.1 
exercise 7.4 
exercise 7.5 
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complex saddle saddle-center 




degenerate saddle real saddle 



Figure 7.4: Stability of a symplectic map in 




generic center degenerate center 



absorb the a factor into x in order to bring tine Henon map for ttie b - -I parameter 
vaiue into the form 

xi+i + x,_i — a — x^j , i — ^, ■■; np , (7.28) 
The 2-dimensional Her\or\ map for b - - I parameter vaiue 

x„+i = fl-x^-y„ 

yn+i - x„ . (7.29) 

is Hamiitonian (sympiectic) in the sense that it preserves area in the [x,y] plane. 

For definitiveness, in numericai caicuiations in examples to foiiow we shaii fix 
(arbitrariiy) the stretching parameter vaiue toa-6, a vaiue iarge enough to guarantee 
that aii roots ofO- f"{x) - x (periodic points) are reai. exercise 8.6 



Example 7.6 2-dimensional symplectic maps: In the 2-dimensional case the 
eigenvalues (5.8) depend only on trM' 

Ai,2 = ^ (trM' + V(trM'-2)(trM' +2)) . (7.30) 

Greene's residue criterion states that the orbit is (i) elliptic if the stability residue |tr M'\ - 
2 < 0, with complex eigenvalues Ai = e'*, A2 = A* = e"'*. If \trM'\ - 2 > 0, A is real, 
and the trajectory is either 

(ii) liyperbolic Ai=e'", A2 = e"'",or (7.31) 
(III) inverse liyperbolic Ai = -e'", A2 = -e"'" . (7.32) 



Example 7.7 Standard map. Given a smooth function g(x), the map 

yn+i = yn+ 8(X„) (7.33) 
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Figure 7.5: Phase portrait for the standard map 
for (a) k = 0: symbols denote periodic orbits, full 
lines represent quasiperiodic orbits, (b) k = 0.3, 
k = 0.85 and k = 1.4: each plot consists of 20 
random initial conditions, each iterated 400 times. 




is an area-preserving map. Tine corresponding ntii iterate Jacobian matrix (4.21) is 

«"(«,«) = n('%;f i). (7.34, 

k-n ' ' 

The map preserves areas, det M = 1, and one can easiiy ciieck tiiat M is sympiectic. 
In particular, one can consider x on tlie unit circle, and y as the conjugate angular 
momentum, with a function g periodic with period 1 . The phase space of the map is 
thus the cylinder Si x R C5 1 stands for the 1 -torus, which is fancy way to say "circle"): 
by taking (7.33) mod 1 the map can be reduced on the 2-torus S2. 

The standard map corresponds to the choice g(x) - kllnsmilnx). When k -0, 
yn+i - yn - yo, so that angular momentum is conserved, and the angle x rotates with 
uniform velocity 



Xn+i = x„ +yo = Jco + (« + l)3'o mod 1 . 



The choice ofyo determines the nature of the motion (in the sense of sect. 2.1.1): for 
yo = we have that every point on theyo - line is stationary, foryo - p/q the motion 
is periodic, and for irrational yo any choice of xo leads to a quasiperiodic motion (see 
figure 7.5(a)). 

Despite the simple structure of the standard map, a complete description of its 
dynamics for arbitrary values of the nonlinear parameter k is fairly complex: this can 
be appreciated by looking at phase portraits of the map for different k values: when 
k is very small the phase space looks very much like a slightly distorted version of 
figure 7.5 (a), while, when k is sufficiently large, single trajectories wander erratically on 
a large fraction of the phase space, as in figure 7.5 (b). 

This gives a glimpse of the typical scenario of transition to chaos for Hamilto- 
nian systems. 

Note that the map (7.33) provides a stroboscopic view of the flow generated by 
a (time-dependent) Hamiltonian 

H{x,y-t) ^\y^+G{x)6y{t) (7.35) 
where 6\ denotes the periodic delta function 



00 

5i{t) = Yj ^(f-'") (7-36) 

/7)=-O0 

and 

G'ix) = -g{x) . (7.37) 

Important features of this map, including transition to global chaos (destruction 
of the last invariant torus), may be tackled by detailed investigation of the stability of 
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periodic orbits. A famiiy of periodic orbits of period Q aiready present in tlie k - Q 
rotation maps can be iabeied by its winding number P/g The Greene residue describes 
tlie stabiiity of a P/Q-cycie: 

RpiQ = \ {2-IxMpiq) . (7.38) 

If RpiQ e (0, 1) tine orbit is eiiiptic, for Rp/g > 1 fhe orbit is liyperboiic orbits, and for 
RpiQ < inverse liyperboiic. 

For k - aii points on tlie yo - P/Q line are periodic with period Q, winding 
number P/Q and marginal stability Rp/q = 0. As soon as k > 0, only a 2Q of such 
orbits survive, according to Poincare-Birkhoff theorem: half of them elliptic, and half 
hyperbolic. If we further vary k in such a way that the residue of the elliptic Q-cycle 
goes through 1, a bifurcation takes place, and two or more periodic orbits of higher 
period are generated. 



7.5 Poincare invariants 



Let C be a region in phase space and y(0) its volume. Denoting the flow of the 
Hamiltonian system by /'(x), the volume of C after a time t is V{t) = /'(C), and 
using (7.24) we derive the Liouville theorem: 



V{t) 



df(x') 



- I dx = \ 

Jf'(C) Jc 

J det {J)dx = ^ dx = 



dx' 
V{0), 



(7.39) 



Hamiltonian flows preserve phase-space volumes. 

The symplectic structure of Hamilton's equations buys us much more than 
the 'incompressibility,' or the phase-space volume conservation. Consider the 
symplectic product of two infinitesimal vectors 



{6x\6x) = 5x^ u)6x = SpiSqi - SqiSpt 



^ {oriented area in the (qi, pi) plane) . 



(7.40) 



i=\ 



Time t later we have 



{6x'\6x) = 6x^ ojJ6x = 6x^a>6x. 



This has the following geometrical meaning. Imagine that there is a reference 
phase-space point. Take two other points infinitesimally close, with the vectors 6x 
and 6x describing their displacements relative to the reference point. Under the 
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dynamics, the three points are mapped to three new points which are still infinites- 
imally close to one another. The meaning of the above expression is that the area 
of the parallelepiped spanned by the three final points is the same as that spanned 
by the initial points. The integral (Stokes theorem) version of this infinitesimal 
area invariance states that for Hamiltonian flows the sum of D oriented areas "V; 
bounded by D loops Q'V,-, one per each {qi,pi) plane, is conserved: 



One can show that also the 4, 6, • • • , 2D phase-space volumes are preserved. The 
phase space is 2D-dimensional, but as there are D coordinate combinations con- 
served by the flow, morally a Hamiltonian flow is D-dimensional. Hence for 
Hamiltonian flows the key notion of dimensionality is D, the number of the de- 
grees of freedom (dof), rather than the phase-space dimensionality d = 2D. 

Dream student Henriette Roux: "Would it kill you to draw some pictures here?" 
A: "Be my guest." 



Resume 

Physicists do Lagrangians and Hamiltonians. Many know of no world other 
than the perfect world of quantum mechanics and quantum field theory in which 
the energy and much else is conserved. From the dynamical point of view, a 
Hamiltonian flow is just a flow, but a flow with a symmetry: the stability matrix 
Aij = oDik H^kj{x) of a Hamiltonian flow Xj = a)ijH j{x) satisfies A^o) + ojA = 0. Its 
integral along the trajectory, the linearization of the flow / that we call the 'Jaco- 
bian matrix,' is symplectic, and a Hamiltonian flow is thus a canonical transforma- 
tion in the sense that the Hamiltonian time evolution x' = f{x) is a transformation 
whose Unearization (Jacobian matrix) / - dx' jdx preserves the symplectic form, 
7^6t»7 = 6t» . This implies that A are in the symplectic algebra sp{2D), and that the 
2D-dimensional Hamiltonian phase-space flow preserves D oriented infinitesimal 
volumes, or Poincaie invariants. The Liouville phase-space volume conservation 
is one consequence of this invariance. 

While symplectic invariance enforces |A| - 1 for complex eigenvalue pairs 
and precludes existence of attracting equilibria and limit cycles typical of dissipa- 
tive flows, for hyperbolic equilibria and periodic orbits |A| > 1, and the pairing 
requirement only enforces a particular value on the 1/A contracting direction. 
Hence the description of chaotic dynamics as a sequence of saddle visitations is 
the same for the Hamiltonian and dissipative systems. You might find symplec- 
ticity beautiful. Once you understand that every time you have a symmetry, you chapter 1 
should use it, you might curse the day [7.26] you learned to say 'symplectic'. 




(7.41) 




in depth: 

appendix C.4, p. 818 
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Commentary 

In theory there is no difference between theory and prac- 
tice. In practice there is. 
— Yogi Berra 

Remark 7.1 Hamiltonian dynamics, sources. If you are reading this book, in theory 
you already know everything that is in this chapter. In practice you do not. Try this: 
Put your right hand on your heart and say: "I understand why nature prefers symplectic 
geometry." Honest? 

Where does the skew-symmetric a) come from? Newton / = ma law for a motion in 
a potential is mq - -dV . Rewrite this as a pair of first order ODEs, q - p/m, p = -dV , 
define the total energy H{q, p) - p^/2m + V(q) , and voila, the equation of motion take on 
the symplectic form (7.3). What makes this important is the fact that the evolution in time 
(and more generally any canonical transformation) preserves this symplectic structure, as 
shown in sect. 7.3.1. Another way to put it: a gradient flow i = -dV{x) contracts a state 
space volume into a fixed point. When that happens, V{x) is a 'Lyapunov function', and 
the equilibrium ,y = is 'Lyapunov asymptotically stable'. In contrast, the '-' sign in the 
symplectic action on {q, p) coordinates, p - -dV induces a rotation, and conservation of 
phase-space areas: for a symplectic flow there can be no volume contraction. 

Out there there are centuries of accumulated literature on Hamilton, Lagrange, Ja- 
cobi etc. formulation of mechanics, some of it excellent. In context of what we will 
need here, we make a very subjective recommendation-we enjoyed reading Percival and 
Richards [7.3] and Ozorio de Almeida [7.4]. Exposition of sect. 7.2 follows Dragt [7.15]. 
There are two conventions in literature for what the integer argument of Sp(- ■ ■) stands 
for: either Sp(D) or Sp((i) (used, for example, in refs. [7.15, 7.17]), where D - dof, and 
d - 2D. As explained in Chapter 13 of ref. [7.17], symplectic groups are the 'nega- 
tive dimensional,' d ^ -d sisters of the orthogonal groups, so only the second notation 
makes sense in the grander scheme of things. Mathematicians can even make sense of the 
d =odd-dimensional case, see Proctor [7.18, 7.19], by dropping the requirement that tj is 
non-degenerate, and defining a symplectic group Sp(A^, w) acting on a vector space M. as 
a subgroup of Gl{M) which preserves a skew-symmetric bilinear form oj of maximal pos- 
sible rank. The odd symplectic groups Sp(2D +1) are not semisimple. If you care about 
group theory for its own sake (the dynamical systems symmetry reduction techniques of 
chapter 10 are still too primitive to be applicable to Quantum Field Theory), chapter 14 
of ref. [7.17] is fun, too. 

Referring to the Sp((i) Lie algebra elements as 'Hamiltonian matrices' as one some- 
times does [7.15, 7.20] conflicts with what is meant by a 'Hamiltonian matrix' in quantum 
mechanics: the quantum Hamiltonian sandwiched between vectors taken from any com- 
plete set of quantum states. We are not sure where this name comes from; Dragt cites 
refs. [7.21, 7.22], and chapter 17 of his own book in progress [7.16]. Fulton and Har- 
ris [7.21] use it. Certainly Van Loan [7.23] uses in 1981, and Taussky in 1972. Might go 
all the way back to Sylvester? 

Dream student Henriette Roux wants to know: "Dynamics equals a Hamiltonian plus a 
bracket. Why don't you just say it?" A: "It is true that in the tunnel vision of atomic 
mechanicians the world is Hamiltonian. But it is much more wondrous than that. This 
chapter starts with Newton 1687: force equals acceleration, and we always replace a 
higher order time derivative with a set of first order equations. If there are constraints, or 



newton - 18jim2012 



ChaosBook.org versionl4.4.I, Apr 21 2013 



CHAPTER 7. HAMILTONIAN DYNAMICS 



138 



fully relativistic Quantum Field Theory is your thing, the tool of choice is to recast New- 
ton equations as a Lagrangian 1788 variational principle. If you still live in material but 
non-relativistic world and have not gotten beyond Heisenberg 1925, you will find Hamil- 
ton's 1827 principal function handy. The question is not whether the world is Hamiltonian 
- it is not - but why it is so often profitably formulated this way. For Maupertuis 1744 vari- 
ational principle was a proof of God's existence; for Lagrange who made it mathematics, 
it was just a trick. Our sect. 33.1.1 "Semiclassical evolution" is an attempt to get inside 17 
year old Hamilton's head, but it is quite certain that he did not get to it the way we think 
about it today. He got to the 'Hamiltonian' by studying optics, where the symplectic struc- 
ture emerges as the leading WKB approximation to wave optics; higher order corrections 
destroy it again. In dynamical systems theory, the densities of trajectories are transported 
by Liouville evolution operators, as explained here in sect. 16.6. Evolution in time is a 
one-parameter Lie group, and Lie groups act on functions infinitesimally by derivatives. 
If the evolution preserves additional symmetries, these derivatives have to respect them, 
and so 'brackets' emerge as a statement of symplectic invariance of the flow. Dynamics 
with a symplectic structure are just a special case of how dynamics moves densities of 
trajectories around. Newton is deep, Poisson brackets are technology and thus they ap- 
pear naturally only by the time we get to chapter 16. Any narrative is of necessity linear, 
and putting Poisson ahead of Newton [7.1] would be a disservice to you, the student. But 
if you insist: Dragt and Habib [7.24, 7.15] offer a concise discussion of symplectic Lie 
operators and their relation to Poisson brackets. " 

Remark 7.2 Symplectic. The term symplectic -Greek for twining or plaiting together- 
was introduced into mathematics by Hermann Weyl. 'Canonical' lineage is church- 
doctrinal: Greek 'kanon,' referring to a reed used for measurement, came to mean in 
Latin a rule or a standard. 

Remark 7.3 The sign convention of oj. The overall sign of a>, the symplectic invariant 
in (7.3), is set by the convention that the Hamilton's principal function (for energy con- 

Pidq, - Et. With this sign convention the action 
along a classical path is minimal, and the kinetic energy of a free particle is positive. Any 
finite-dimensional symplectic vector space has a Darboux basis such that u) takes form 
(7.9). Dragt [7.15] convention for phase-space variables is as in (7.2). He calls the dy- 
namical trajectory jcq ^ x{x(i, t) the 'transfer map,' something that we will avoid here, as 
it conflicts with the well established use of 'transfer matrices' in statistical mechanics. 

Remark 7.4 Loxodromic quartets. For symplectic flows, real eigenvalues always 
come paired as A, 1/A, and complex eigenvalues come either in A, A* pairs, |A| = 1, or 
A, 1/A, A*, 1/A* loxodromic quartets. As most maps studied in introductory nonlinear 
dynamics are 2d, you have perhaps never seen a loxodromic quartet. How likely are we to 
run into such things in higher dimensions? According to a very extensive study of periodic 
orbits of a driven billiard with a four dimensional phase space, carried in ref. [7.28], the 
three kinds of eigenvalues occur with about the same likelihood. 

Remark 7.5 Standard map. Standard maps model free rotors under the influence of 
short periodic pulses, as can be physically implemented, for instance, by pulsed optical 
lattices in cold atoms physics. On the theoretical side, standard maps illustrate a number 
of important features: small k values provide an example of KAM perturbative regime 
(see ref. [7. 11]), while larger A:'s illustrate deterministic chaotic transport [7.9, 7.10], and 
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the transition to global chaos presents remarkable universality features [7.5, 7.12, 7.7]. 
The quantum counterpart of this model has been widely investigated, as the first example 
where phenomena like quantum dynamical localization have been observed [7.13]. Sta- 
bility residue was introduced by Greene [7.12]. For some hands-on experience of the 
standard map, download Meiss simulation code [7. 14]. 
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Exercises 



7.1. Complex nonlinear Schrodinger equation. Con- 
sider the complex nonlinear Schrodinger equation in one 
spatial dimension [7.26]: 



deb d^d) , 



(a) 









1 














1 


-1 








-c 





-1 


c 






Show that the function i/' : R ^ C defining the 
traveling wave solution (^(.x, t) - il/{x-ct) for c > 
satisfies a second-order complex diflTerential equa- 
tion equivalent to a Hamiltonian system in rel- 
ative to the noncanonical symplectic form whose 
matrix is given by 



Wr — 



(b) Analyze the equilibria of the resulting Ha- 
miltonian system in R'* and determine their linear 
stability properties. 

(c) Let = e"''^'^fl(i) for a real function a{s) and 
determine a second order equation for a{s). Show 
that the resulting equation is Hamiltonian and has 
heteroclinic orbits for /3 <Q. Find them. 

(d) Find 'soli ton' solutions for the complex nonlinear 
Schrodinger equation. 

(Luz V. Vela-Arevalo) 

7.2. Symplectic vs. Hamiltonian matrices. In the 

language of group theory, symplectic matrices form the 
symplectic Lie group Sp(ii), while the Hamiltonian ma- 
trices form the symplectic Lie algebra sp{d), or the al- 
gebra of generators of infinitesimal symplectic transfor- 
mations. This exercise illustrates the relation between 
the two: 

(a) Show that if a constant matrix A satisfy the 
Hamiltonian matrix condition (7.14), then J{t) - 
exp(fA) , f e R, satisfies the symplectic condition 
(7.9), i.e., 7(f) is a symplectic matrix. 

(b) Show that if matrices To satisfy the Hamiltonian 
matrix condition (7.14), then g((p) = exp(0 ■ T) , 

i.e.. 



, satisfies the symplectic condition (7.9), 
I is a symplectic matrix. 

(A few hints: (i) expand exp(A) , A - cp - T , as a power 
series in A. Or, (ii) use the linearized evolution equation 
(7.18).) 

7.3. When is a linear transformation canonical? 



(a) Let A be a [n X n] invertible matrix. Show that 
the map (p : R^" R^" given by (q,p) 
(Aq, (A ')^p) is a canonical transformation. 

(b) If R is a rotation in R^, show that the map (q, p) i-» 
(R q, R p) is a canonical transformation. 

(Luz V. Vela-Arevalo) 

7.4. Determinants of symplectic matrices. Show that 
the determinant of a symplectic matrix is H- 1 , by going 
through the following steps: 

(a) use (7.27) to prove that for eigenvalue pairs each 
member has the same multiplicity (the same holds 
for quartet members), 

(b) prove that the joint multiplicity of /I = + 1 is even, 

(c) show that the multiplicities of A - I and A - -I 
cannot be both odd. (Hint: write 

P(A) = (A-lf'"*\A+ if'^'^QiA) 

and show that g(l) = 0). 

7.5. Cherry's example. What follows refs. [7.25, 7.27] is 
mostly a reading exercise, about a Hamiltonian system 
that is linearly stable but nonlinearly unstable. Consider 
the Hamiltonian system on R'* given by 

H = ^(qi + Pi) - {qj + pI) + ]^P2{p\ - q\) - qiqipi- 

(a) Show that this system has an equilibrium at the 
origin, which is linearly stable. (The linearized 
system consists of two uncoupled oscillators with 
frequencies in ratios 2:1). 

(b) Convince yourself that the following is a family of 
solutions parameterize by a constant r: 



r-COS(f-T) 

q\ = - V2— , 

t - T 

-sin(f - r) 



t - T 



qi = 



pi 



cos 2(f - t) 

f-T 

sin 2(f - t) 

f-T 



These solutions clearly blow up in finite time; 
however they start at f = at a distance V3/t 
from the origin, so by choosing t large, we can 
find solutions starting arbitrarily close to the ori- 
gin, yet going to infinity in a finite time, so the 
origin is nonlinearly unstable. 

(Luz V. Vela-Arevalo) 
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